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r-| , Supersymmetric Grand Unified Theories have a rich spectrum of particles barely heav- 
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ier than the intermediate vector bosons. As their non-supersymmetric counterparts, they 
lead to many relations among low energy observables. But the precise form of the pre- 
dictions is modified by the extended spectrum. If the masses of these new particles are 
comparable to Mz, the standard computation of their effect becomes inaccurate. We 
present a detail discussion of the correct procedure, and carry out the relevant computa- 
tions to one loop order. The procedure we propose has the advantage over other existing 
methods that two-loop running of gauge couplings can be incorporated readily and con- 
sistently. Attention is paid to the special treatment that the top and Higgs particles must 
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receive. The size of the effect is explored for a range of parameters in the minimal super- 
symmetric SU (5) grand-unified theory with radiative breaking. It is found that the naive 
(leading-log) computation can be fairly inaccurate. 
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1. Introduction 

Grand Unified Tlieories possess many attractive features. They simplify enormously 
the description of matter content of the standard model by reducing the number of irre- 
ducible representations of the gauge group needed to account for all observed particles. 
They explain electric charge quantization. They bring all non-gravitational interactions 
under a common umbrella treatment. Yet, the minimal standard model extension to an 
SU(5) GUT is all but ruled out: not only does proton not decay at the predicted rate, 
precise measurements of the gauge couplings of strong and electroweak interactions have 
led to the observation that at no scale do they unify. 

This rather grim observation is, nevertheless, not generic of GUT theories, and it 
seems premature to dismiss the class of theories on the basis of the failure of the minimal 
version. Remarkably, the minimal supersymmetric extension of the minimal SU(5) GUT is 
free of the aforementioned problems. That the three gauge couplings unify |0] in that case is 
not a surprise: an additional parameter, Msusy, is introduced which can be chosen so that 
couplings do unify. What is remarkable is that unification occurs at a scale Mqut high 
enough that proton stability is not in confiict with observation, yet bellow the Planck scale, 
making it plausible that calculations with presently understood techniques are sensible. 
Moreover, the additional parameter, Msusy, which describes the scale at which the running 
of the gauge couplings changes from what is dictated by the standard model to what is 
dictated by its supersymmetric extension, comes out rather small, just about equal to the 
scale of electroweak symmetry breaking. If one interprets MgusY as some average mass 
of the supersymmetric partners of standard particles, the inescapable conclusion is that a 
rich spectrum of new elementary particles is awaiting discovery in the few hundred GeV 
range (well accessible to planned next generation hadron colliders). 

One ought not to rush into conclusions, though. Barbieri and Hallf^ have pointed 
out that the standard analysis of running of gauge couplings assumes all the superheavy 
particles (those with masses naturally of order of Mqut) are degenerate. The value of 
Msusy? they argue, can be shifted if one drops this assumption and includes the effects of 
the non-degenerate superheavy thresholds. On closer examination, it has been noticed 
that in the minimal SUSY SU(5) model the dependence of Msusy on superheavy thresholds 
is rather weak. 

For example, if as{Mz) = 0.125, the value of Msusy exclusive of superheavy threshold 
effects comes out to be well below the electroweak scale. Superheavy thresholds can push 
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-^SUSY back up. But to obtain MgusY > Mz one must have that the superheavy super- 
multiplets that transform as (8, l)o + (1, 3)o under the gauge group SU (3) x SU (2) x U (1) 
acquire masses some eight orders of magnitude below Mqut- Worse yet, to keep Mqut 
weU below the Planck scale one must further complicate the spectrum of superheavies, 
splitting the masses of these two multiplets by two orders of magnitude. The simple com- 
pelling picture of a dessert between the electroweak and GUT scales is lost. Instead one 
must populate the dessert with particles at two new intermediate scales. For the remaining 
of this work we neglect these heavy threshold effects. 

In addition to superheavy threshold effects, one must consider the effects of light 
thresholds (non-degeneracy of the SUSY particles with masses naturally of order the elec- 
troweak scale). Ross and Roberts^] argue that because in most models the average mass 
of the colored superpartners is larger than that of the uncolored ones, a simple param- 
eterization in terms of a single threshold at MgusY may be misleading. Langacker and 
Polonski0 point out that one need not interpret Msusy as some sort of average mass for 
the spectrum of SUSY particles. The spectrum of the model is not expected to be simple. 
One can imagine accounting for the complicated spectrum in some accurate fashion, and 
then finding a scale Msusy which would mimic the effects of a complicated spectrum. Only 
this scale does not have a clear physical meaning. In fact, they find it possible to construct 
examples of rather heavy SUSY spectra (compared with Mz), yet having Msusy Mz- 

The standard analyses of light threshold effects consists of modifying the evolution 
equations for the gauge couplings at each consecutive threshold]^. More specifically, start- 
ing from unified gauge couplings at Mqut, ai(MGUT) = 0^2 (Mqut) = as (Mqut), one 
uses the renormalization group equation (RGE) to evolve these couplings down to the first 
threshold. The evolution equations are modified by reducing the number of active degrees 
freedom that affect the running, and the couplings are run further down to the next thresh- 
old, and so on. In what follows we shall refer to this method as the 'run-and-match' or 
the 'naive' method. This method is justified, and expected to be a good approximation, 
provided all of the SUSY particles are much heavier than the Z-boson. Also, this kind of 
analyses by necessity assumes unbroken electroweak symmetry, which again is reasonable 
at mass scales well above Mz- 

In this paper we address the issue of how to compute these light threshold effects when 
the SUSY particles are not necessarily much heavier than the Z-boson. The motivation 
for this comes not just from the discussion just given, but from the recent observation that 
phenomenologically viable models that incorporate SUSY SU(5) and radiative electroweak 
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breaking invariably contain many particles with masses barely exceeding Mz- Moreover, 
we will craft the method of computation so that two-loop running effects are incorporated 
consistently. 

We will first formulate carefully the question we want to address. This is done in 
sections 2 and 3 where we give the fundamental equations that allow us to relate the 
measurements of gauge couplings as reported by, say, LEP/SLC experiments to the gauge 
couplings of the SUSY model. The important discussion of the effects of two-loop running 
can be found in section 4. The difference between a calculation of threshold effects based 
on effective field theories versus calculations based on, say, the so called "star scheme" is 
discussed there. In section 5 we develop general formulae for threshold effects modifications 
to relations among parameters dictated by grand unified theories. In section 6 we exhibit 
the result of one loop computations of the quantities that go into the general formulae of 
the preceding section. Special care must be taken in dealing with the Higgs sector and 
the top quark, and this is addressed in section 7. In preparation for the discussion and 
presentation of our numerical analysis of section 9, we review the minimal supergravity 
SU(5) model in section 8. Our conclusions can be found in section 10. 

There is a long history of calculations of radiative corrections in and beyond the 
electroweak theory]^. Calculations in supersymmetric theories date back to refs. ||^. No 
doubt the trivial calculations of one loop diagrams in this paper can be found elsewhere. 
The emphasis in this paper is on the treatment of light threshold effects in grand-unified 
theories, and in particular in SUSY-GUTs. It may not be obvious to the reader but, as we 
will see (sect. 3), the computations involved coincide with those of radiative corrections to 
electroweak parameters. Thus, for example, while ref. presents a detailed discussion of 
radiative corrections to electroweak parameters in SUSY extensions of the standard model 
(see table 1 and figs. 39-45 of that work), the same work treats the case of grand-unified 
SUSY theories as containing a single common mass threshold at a scale '/x' (see figs. 47-49 
there). We believe ours is therefore the most complete treatment of the light thresholds 
in SUSY GUTs to date. Moreover, we have failed to find as detailed a presentation of 
the numerical size of the effects for the minimal SUSY SU{5) GUT theory with radiative 
breaking as we give here (c/, sect. 8). 
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2. Threshold Corrections 



Let us define what is meant by threshold corrections in this paper. The gauge sector 
of the standard model of electroweak and strong interactions has four independent pa- 
rameters: the three gauge couplingsB (71,2,3 and the vacuum expectation value v. These 
parameters are not all independent in a unified theory. For example, in minimal SUSY 
SU(5) GUT one has 

sm 9=- + — , 2.1 

5 15 as 

where sin^ ^ = gf/igf + (/I)? cte.m. = 9i92/^'^{9i + 92) and as = 9s/4:7i are all defined 
in terms of the running coupling constants, and the relation holds as a function of the 
renormalization scale /x. This relation holds in the leading log approximation (corrections 
are suppressed by 1/ log(MGUT//i))- It assumes a spectrum of particles much lighter than 
-^GUT, the "light particles", consisting of the particles in the standard model with two 
higgs doublets, and adding the superpartner of each field. 

To test whether this relation holds one must carefully define the terms that enter into 
it. It would seem appropriate to cast this relation in terms of physical observables. For 
example, one may write 

= 1 + 1^ (2.2) 
5 15 as 

or 

1-M^/M| = 1 + -^^. (2.3) 
5 15 as 

These are not equivalent, and cannot both hold. The discrepancy comes from the differ- 
ence, beyond tree level, between two quantities that agree at tree level. Given a precise 
definition of the four observables G^, Mz, ae.m. and as, the equations can be corrected to 
account for the difference: 





1 - 1 = - + 1^ + 5^ (2.4a) 

^G^Ml) 5 15 as ^ ' 

1 - Ml, /Ml = \ + L?^ + S2 {2 Ah) 
5 15 a c 



^ We assume, for definiteness, that dimensional regularization is used and that lagrangian 
parameters are defined through minimal subtraction. 



5 



Contributions to 81^2 fall into two distinct classes: those that are there in the standard 
model, and those that arise from the virtual effects of new particles beyond the standard 
model's. It is the latter that we would like to focus on and which give rise to what we 
call "light threshold effects". The idea is that one need not recompute the effect of the 
standard model particles every time an analysis like this is performed. Moreover, one needs 
not grope with the issue of which physical observables to use to represent sin^ 9. Instead, 
one may perform a standard model based analysis to extract the values of the couplings 
fi'1,2,3 and V. If there is an underlying SUSY SU(5) GUT, the couplings extracted from an 
analysis based on the full spectrum are then expected to satisfy a relation like in eq. ( |2.1|) . 
If only the standard model's spectrum is used in the extraction of the couplings, there will 
be a small correction 

sin2^=^ + -^^ + flight, (2.5) 
5 15 

where 5iight is precisely what we intend to mean by "light threshold corrections" . Eq. ( |2.5| ) 
is used in ref. to implicitly define light threshold effects. The arguments leading to Jught 
there are different and simpler, given the nature of the approximations in that paper. Here 
we intend to generalize the results of ref. 0, and we will find it useful to check that our 
results for 5iight reduce to those of ref. in the appropriate limit. 

3. Generalities 

Gauge coupling constants are not physical observables. When an experiment reports 
on the value of a gauge coupling, there is an implicit translation of some observables into 
these theoretical constructs. This is not to say that the values of gauge couplings lack 
in importance. On the contrary, given a well specified definition of these couplings, they 
encode concisely the results of measurements. 

When an experimental result is analyzed to extract the value of gauge coupling con- 
stants, the actual value obtained depends in detail on the theoretical assumptions. Given 
two different models with the same gauge group but different particle content, as is the 
case of the standard model and its supersymmetric extension, the extracted values of gauge 
couplings from the same set of observables is a priori different for the two models. In prin- 
ciple one could analyze the observables directly under the different set of assumptions, ze, 
different models, and thus extract the values of gauge coupling constants appropriate to 
those assumptions. In practice however, it is impossible to reanalyze the experiment for 
each new set of assumptions. 
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The results of experiments are therefore most often given in terms of an analysis 
based on the standard model. This makes sense. The standard model is the most concise 
model of elementary interactions consistent with all present observations. Clearly what 
we need is a means for translating experimentally determined values of gauge couplings in 
the standard model into gauge couplings in extensions of the standard model. 

It is possible to have such a translation. An obvious case is one in which the model at 
hand is an extension of the standard model that has only very heavy particles with hard 
masses. By this we mean that their masses are much larger than the energies at which 
the experiments are carried out, and that the masses are symmetric under the electroweak 
SU{2) X U{1), up to small corrections. Then the decoupling theorem guarantees that the 
physics at experimental energies is described by an effective Lagrangian which corresponds 
precisely to the standard model. The effective coupling constants of this effective theory are 
numerically equal to the standard model's. The gauge couplings of the underlying theory 
are easily related to those of the effective theory, by matching conditions at the scale 
of the heavy particles. Non-supersymmetric GUTs furnish an example. The underlying 
coupling constant, aouT, can be derived by matching at Mqut, so one must determine 
experimentally the electroweak and strong couplings at a low scale, and then run them up 
to Mqut where the effective theory can be matched to the full underlying theory. 

In fact, this method can also be used even when the heavy particles are only moder- 
ately heavier than the energy scale of experiments. This is because this method accounts 
for all effects that are logarithmic in the ratio of the experimental scale E and the heavy 
mass M, lii{E/M); corrections are suppressed by powers of E/M, or even of (E/M)'^. 
Therefore it is quite sensible to use this method to analyze SUSY GUTs, especially if one 
is interested in a spectrum of heavy particles in the TeV range. In fact, this is precisely the 
way in which the low threshold effects have been accounted for by many in the past []5|— |7|,pT[] . 
We refer to this method as the 'run-and-match' or 'naive' calculation. 

But one may ask what the proper procedure is, and even whether one exists, for the 
case in which the heavy particles are only barely heavier than the experimental scale, ie, 
E/M is only slightly smaller than unity. After all, many recent analysis of SUSY GUTS 
show that models typically contain several particles in the few hundred GeV range. And 
not only is E/M only slightly smaller than 1; since the multiplicity is high, one expects to 
find large coefficients in front of the order [E/M) correction. 

It turns out it is easy to construct the proper procedure. We will assume that particles 
are heavy enough that they cannot be directly produced at present (else, the theory should 
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be ruled out or confirmed experimentally in short order!). This is important because it 
implies immediately that the only observables available are those of the standard model. 
In other words, we can consider, say, cross sections or decay widths of standard model 
particles for comparison with experiment, and not consider any novel processes of the 
extended theory. Let Oi, i = 1, 2, 3, . . . stand for a collection of observables. For example, 
this can be taken to be the forward-backward asymmetry at the peak of the Z-resonance, 
the differential cross section da{e^e~ — » fx'^ fx~)/d cos 9 at a collection of specific center 
of mass energies and angles, and similar cross sections into other leptons or quarks. The 
standard model gives expressions Fi for these observables in terms of its parameters Qk, 



Now, any extension of the standard model will also give expressions Fi for the same 
observables in term of its own parameters. Among these parameters are those which 
have a direct correspondence to those of the standard model, gk, like gauge couplings and 
particle masses. In addition there are new parameters e„ characterizing the new physics, 
like new particle masses. So one has. 



Although conceptually the same, one must differentiate between the couplings gk and gk- 
In fact, it is this difference that we are trying to calculate for the case of gauge couplings! 
In other words, the problem at hand is to calculate gk from knowledge of gk, given an 
assumption on the values of the parameters e„. 

Our task is to invert eqs (|3.1|) -(|3.2|) for gk in terms of gk (and, implicitly, Cn). Now, 
the additional degrees of freedom enter into the expressions for observables as heavy virtual 
particles. Therefore, their effects are small. We can write 



where 5gk are small one loop order corrections which are expressed as functions of gk (and, 
implicitly, e^). It turns out to be convenient to write the functions F as a sum of two 
pieces, one identical to the standard model's, plus the rest: 







gk = gk + Sgk , 



(3.3) 



Fi{gk, en) = Fi{gk) + AFi{gk, e„) 



(3.4) 
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This is convenient because AF can be regarded as a small quantity, the same order as 5g. 
This means also that we are assuming the same renormalization scheme and gauge choice 
in both models. Combining equations and keeping only leading order in small quantities, 
one obtains 
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/ J9i^—i9k) + AFi{gk,en) = . (3.5) 
^ ogi 



Equation (|3.5| ) is a set of simultaneous equations for 5gk- In the first term, the derivative 
should be evaluated at tree level, because the coefficient 5gi is already of one loop 
order. This is a welcome simplification: we only need to compute the one loop effects of the 
new particles, because the standard ones have already been included in the experimental 
determination of the couplings gk- 

It is important to note that this discussion depends crucially on approximate de- 
coupling. Were the new particles not heavy enough, the dependence of observables on 
kinematic variables (like s and t in cross sections) would differ from one model to the 
other. To the extent that deviations from the standard model cannot be inferred from 
kinematic dependence of observables one can neglect these kinematic effects in AF. Al- 
though this observation seems rather innocuous, it is of practical importance. We need it 
to ensure compatibility of expressions ( p.l|) and (|3.2|) . 



Now, it would seem the program of determining the threshold corrections 6gk is com- 
plete: calculate the functions Fi at tree level, and the virtual effects of the new particles 
at one loop AF^, and plug into eq. ( p.5| ). Then solve for dgk- Still, this program involves a 
choice of specific observables, labeled by i. Recall that we would like to be able to connect 
the (jk to the gk without detailed knowledge of the experimental data. 

The Feynman diagrams that contribute to the functions AFi can be classed in two 
groups: corrections to standard particle propagators and vertex or box corrections. The 
former are process independent, ie, they enter in a universal form into the relation between 
g and g. The latter depend on each process, but are often smaller than the former. We 
therefore drop them. In standard lingo, we retain only oblique corrections. In any case 
this kind of truncation is implicitly done in the run-and-match approach. Our analysis 
is intended as an improvement on that approach, but it is clear that it is incomplete in 
this regard. (It should also be clear, from the above discussion, that short of a complete 
reanalysis of experimental data, this is the best one can do). 
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4. Two Loop Running 

The relation between coupling constants in eq. ( |2.1| 



5 15 as 

follows from the solution to the renormalization group equations at one-loop order 



namely 



^ ^ -6oln4- (4.2) 



At two-loops the renormalization group equations are difficult to solve. They are a set of 
coupled equations. Neglecting yukawa couplings one has 

^^=^f^^7f + E^S"'^'^|- (4-3) 

j 

How does this modify the relation between gauge couplings in eq. (|2.1D ? In principle there 
is an analytic expression that expresses the new relation between the couplings. In practice 
one recognizes that, because the gauge couplings are small, the corrections to the relation 
in eq. ( p.l| ) are correspondingly small. So one may write 

sin^ ^ = 1 + 1 ^ + 5two-ioop . (4.4) 
5 15 tts 

It is customary to calculate the correction 5two-ioop by solving the renormalization groups 
eqs. ( ^.31 ) numerically and comparing to the solution of the one- loop version, eq. ( [4.2| ) . 
There are two important points that must be stressed: 

(i) The combined correction from light threshold and two-loop running effects is computed 
by simply adding together the corresponding corrections, 5iight and 5two-ioop; and 

(ii) Both corrections (5iight and 5two-ioop) are formally of the same order, namely 0{ai). 
The first of these points is not only obvious but holds to arbitrary order. The low 

energy effective theory has gauge couplings which satisfy renormalization group equations. 
In principle, these may be solved to arbitrary accuracy. In the effective theory one can 
then make predictions of physical observables, and carry out the computations to the same 
accuracy as the running. Effective field theories are particularly good at making this sort 
of logic, and the corresponding calculation, transparent. 
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Consider, in contrast, the prediction using the "star" scheme [1T0| . One is instructed 
to compute star couphng constants gi*{q'^) in terms of the bare couphng constants Qi bare, 

7T^-72^ = n.(^^) (4.5) 

yi*iy ) iJi bare 

where the Iii{q^) are appropriately chosen two-point functions, directly related to self- 
energies of gauge-particles. It is easy to see how the star and effective field theory ap- 
proaches give the same answers at one loop. For example, the right sides of eq. (|4.5| ) is 
dominated by a large logarithm when is very small compared to the cut-off, A^, and 
this logarithm precisely corresponds to the one in eq. ( |4.2D if we make the correspondence 

\q^/K^\^nV^^l 

In the star scheme a two loop contribution to Ili{q^) is necessarily of order g'f. There- 
fore, it gives a correction Sgf/gf of order gfg'j. This is in contrast with the correction 
(^two-loop, which gives 5gf / gf ~ g'j- To see this, we may solve a simpler equation that 
does the power counting correctly. Take, for example, pure QCD, which corresponds to 
neglecting the couplings gi^2 in the running of (73 in eq. (|4.3|) , 

^ ^7 = ^^9l + h,gl . (4.6) 
The solution is well known, and has a simple expression if one retains only leading terms 
in 53, 

1 1 A 1 ^0 , ^1 1 f3(^0)^ . . „^ 

- 60 In ^ + 7- In — —-2- . (4.7) 



The second term on the right side cannot result from the two loop computation of the 
two-point function 11 in eq. ( [4.5|) . It is not clear to us how the star scheme can incorporate 
the two-loop running, if at all. 

That flight and 5two-ioop are of the same order is readily seen by explicit computation. 
Indeed, an expression for ^Hght appears in the next section, and is explicitly of order a. 
And comparing eqs. ([4.2|) and (f4.7| ) one infers the correction to the running coupling 6g'^, 
with 5(l/fi'^) = -Sg'^/g'^ ~ 0{g'^), or Sg"^ ~ g"^, which leads to 5two-ioop ~ «. 

It is worth noting that formally the heavy-threshold effects 5heavy are of the same 
order as ^ught and 5two-ioop. i] 

Numerical analysis of the effects of two loop running (including yukawa couplings) has 
been extensive |]12|. Our numerical results do not include the effects of two-loop running. 
Our emphasis is on the light threshold effects and how they depend on the spectrum. The 
effects of two-loop running are independent of the details of the spectrum. The interested 
reader may incorporate the two-loop running effect simply by adding an overall constant 
(ie, spectrum independent), ^two-loop, to our predictions for sin^ 6*. 
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5. Derivation of Master Formulas. 



We are ready to proceed with explicit calculations. Although the expression for the 
threshold corrections Sgk will be process independent, it is necessary to choose particular 
(gedanken) process to obtain it. 

We first consider a particularly simple example, /uz/^ ez/e scattering. We hope this 
will serve to illustrate the general approach with a minimum of unnecessary complications. 
We take the observable O as the cross section (at some kinematic point), divided by the 
phase space, ie, the square-modulus of the invariant amplitude, 

= F{g)= 1^1' (5.1) 

helicities 

At tree level this is given by the Feynman diagram in fig. |l| 

Now, we really only need to extract the dependence on couplings, without paying 
careful attention to the kinematic factors. Therefore we write 

helicities z / helicities 

Here s is the standard kinematic variable, g2 is the SU{2) coupling constant and v is the 
electroweak breaking vacuum expectation value. Clearly A is just a product of spinors and 
gamma matrices. 

The one loop virtual effects of new heavy particles includes vertex, box and propagator 
corrections; these are shown in fig. ||. As discussed previously, we will neglect the vertex 
corrections. The propagator corrections are conveniently expressed in terms of the self- 
energy (defined as the IPI amputated two point function for the charged vector boson): 

^Wfiuiq) = {q^9^,v - ?/xQ^)n{J^ + q^,qu^w^ (5.3) 

In terms of these the VF-boson propagator in Feynman-'tHooft gauge is 

. 9^^u-q|^q./q^ . q^^qu/q^ 



The function AF that we need to compute is defined by eq. ( p.4|) , which tells us that 
we should take the difference between in the extended model and in the standard 
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model. The tree level amplitudes are equal in both models, so the difference comes from 
the cross product of the tree level amplitude and the one-loop amplitude: 

" (T) 

AF^29.^Re ^"'^ J] \A\K (5.5) 

^ ^2 ) helicities 

Here by we mean that part of the self-energy which is not already in the standard 
model calculation, ie, which comes from the new particles of the extended model running 
around in the loop. Using eqs (p.l| ) - (|5.5|) in eq. ( ^75|) , one obtains 



r 9 ^ d 

og2T. '^(>v — 

og2 ov 



^'w) + TT^^^^^^^w = ^ ■ (5.6) 



or, simplifying, 

^ + ^^^ = -RenmW (6,7) 

Similar expressions can be obtained for other couplings by considering other processes. 
From neutrino scattering one obtains 

^^^^-HefirO). (5.8) 

while from quark scattering 

5gl = -glReIif\Ml) . (5.9) 

Herei gi and g^, stand for the gauge couplings for the U{1) and SU (3) groups, respectively, 
and liz and lig are the Z-boson and gluon self energies. Again, the hat over these quantities 
reminds us to include only non-standard model contributions. 

At first sight these expressions appear nonsensical: the dependence on the kinematic 
variable, s, is manifestly different on both sides of the equation. Were high precision 
experiments performed over a wide range of energies, one could either verify or rule out 
these effects. In practice precision experiments are available only at very low energies, 
s ^ M^, or on resonance, s ~ or M^. We will not assume that the functional form 
of (s) or (s) can be well approximated by truncating their expansion to contain 
only a pole and constant terms. Although this would match the functional form on both 



^ In this section, the normalization of g\ is the standard model one. In the following sections 
we will switch, without warning, to the normalization that is more appropriate for GUT theories. 



That is, we will absorb a factor of a/5/3 in g\. 
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sides of eqs. (|5.7| ) and (|5.8|) , one of our objectives is to retain the full effects of these 
self-energies. Thus, from eqs. ( p.7|) and ( p78|) one obtains a set of four equations for three 
unknowns: 



92 ^ 

2 

+ ^ = -ReUPiMi:) (5.106) 



{9f + 9l) 



dv s 



r 2 

^ = -hm-^Renl.^Vs) (S.lOd) 

Clearly, the last two equations are not generally mutually compatible. It is easy to see 
that the difference between the right hand side of eqs. ( ^.10(^) and ( ^.10d|) is a contribution 
too Veltman's p parameter. There are two ways of dealing with this issue. Since there is no 
direct evidence for deviations from the standard model, one may consider such contribution 
to p negligible. This is in keeping with the general philosophy of our method. Nevertheless, 
one may object that the experimental determination of the p parameter is not at the 
same level of precision as, say, measurements of Mz or forward-backward assymmetries. 
Therefore, our theoretical results for, say, sin^ 6 would inherently carry an error which 
could be much larger than the precision with which it is measured. Alternatively, one 
may introduce the p parameter itself into the analysis. This corresponds to replacing 
eqs. (fUJBj ) and (fUJ^ by 

f±f .^.^^-Hen-(Ml) (5.10.) 

^ + ^ = -\iinA.ReUP(s) (5.10d') 

Although p is not a fundamental parameter of the standard model, it is usually included 
in experimental fits. 

We will investigate both approaches. One cannot choose one over the other as a matter 
of principle: it is the way that the experimental fit is performed that determines which 
analysis should be chosen. We present both for comparison. Clearly it is of interest to 
carry out the same type of analysis with the extended set of six parameters of Golden and 



Randall ||19|| and Altarelli and Barbierip^ — adding not just p to the standard model's. 
We will return to such analysis in a forthcoming publication. 
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In what follows we will refer to the approaches that neglect and include the contribu- 
tions to p as method-I and method-II, respectively. 

We are now ready to compute threshold corrections 5iight as defined in section 2, 
eq. ( ^.51 ). Using the results above, one obtains in method-I 



^light 



dp 



207r 



-1 + 2sin^^ - 5sin^^ 



5 sin^ e -I 



1 



sm i) sm i) o 



(5.11) 



where we have introduced the reduced self-energies 
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Renf)(M|) 



Renf)(M|) 
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e.m. 
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47r sin^^^ 

Cle.m. 1 



47r sin^ 



47r sin^ 6'cos2 6'^^ 



(5.12a) 
(5.126) 

(5.12c) 
(5.12d) 



These are functions of the masses of the virtual particles. When any one such mass becomes 
large, the leading contribution to these functions is in the form of a logarithm of the ratio of 
the mass to the Mz or Mw- Thus one recovers the form of the result of the usual method, 
described earlier, of matching couplings and changing evolution equations at thresholds. 
We have checked that our expressions agree with the standard run-and-match formulas 
when we take the large mass, SU{2) symmetric limit. 
The threshold corrections in method-II are 



"light 
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-l-f2sin2^-5sin^^ 
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$w - + —^i'^z - $z) + 



sin*^ 



(5.13) 



where in addition to the definitions in (|5.12| ), we have introduced 



lim sReftP' 



47r sin 6 cos^ 6 



(5.14) 



6. Explicit Computations 

In the previous section we set up the stage for explicit computations. In this section 
we present explicit results of one-loop computations. 
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Lets consider first the case of scalar particles contributions to The two required 
one-loop Feynman diagrams are shown in fig. ^. The calculation is performed using di- 
mensional regularization and minimal subtraction. The result can be written succinctly in 
terms of the function 

/i TOi,m2,Q = ^In^ + ^In^ - - -2 / dx ^ In 6.1 

[i^ 3 Jo 1 A* 

where 

M'^{x) = m\x -f m^(l — x) — (^x{\ — x) 

and \x is the renormalization point. For a colored complex scalar of mass m in the R 
representation of color S'[/(3), with Casimir invariant C{R), one has 

^g = C{R)Ii{m,m,Mz) . (6.2) 

The all important effect of mixing enters the computation of and $VK- Writing 
the scalar field component of the neutral current as 

9ij <Pid^(pj , (6.3) 

and the corresponding scalar masses as rriij, one obtains the Z self-energy 

$z = bzjf /i(to^,to„M^) . (6.4) 

The coupling gij includes the usual group theoretic factors and gauge couplings, and the 
angles from the mixing matrix. Consider, for a definite example, a scalar of mass m which 
does not mix, that is, for which the mass and weak eigenstates agree. We only need 
consider scalars in representations which are singlets or doublets under weak-SU{2). Let 
stand for the third component of weak-isospin and Q for electric charge. One obtains, 

$z = (t^ - sin^ eQfli{m, m, Mz) . (6.5) 

The expression for $VK is entirely analogous to that of $z in eqn. (|6.3|). Lets consider the 
case of no mixing. = for weak-singlets, while if the members of the doublet have 
masses mi and m2 

= \h{mi,m2,Mw) . (6.6) 
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Also, 



nil + ^2 + 



2 2m^m2 



In 



(6.7) 



Next we turn to the case of fermions. There is only one Feynman diagram contribution 
to at one-loop; see fig. ^. In the case of fermions, three functions are needed to describe 
the results: 



/2(mi, 7712, q) = 8 / (ix x(l — x) InM^(x)/^^ (6.8a) 
Jo 

13(1711,7712, q) = — ^ / dx [rrioX + m1{l — x)] In M'^{x) / jj,'^ (6.86) 
Q Jo 



T I \ 4mim2 , , „^2/ N / 2 
l^[mx,m2,q) = 5 — / ax InM [x)l\i 



r Jo 

In terms of these, one has for colored fermions 



(6.8c) 



^g = C{R)l2{m,7n,Mz) . 



(6.9) 



We will be generally concerned with right handed, as well as left handed, weak-doublets. 
Therefore it is useful to introduce some additional notation, allowing for the possibility of 
mixing in the neutral current: 



2 sin 6 cos 6 



(6.10) 



These couplings, and g^, involve both mixing and the usual — sin^ 9Q factors. In 
terms of these, 

= i(|4P + |(7gP)(/2(m,,m,,Mz) +/3(m,,mj,Mz)) +Re(4(7j?*)/4(m,,m,,Mz) 

(6.11) 

The expressions for and $^ can in fact be read off from this, but we give them 
explicitly for future reference. For a fermion doublet with components of mass mi and 

= liMmi, 1712, Mw) + himi, 7712, Mw)) (6.12) 
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and 

r 2 2 1 n 

(6.13) 



q 

= li™; 7772" (^2(^1, m2, q) + h{mi, m2, q)) 



2 2 1 

4 1 ^2 4 1 ^1 , / 4 4 
In — r — mi m — r- H — (mi — m. 



We have neglected mixing, which is a good approximation in the quark sector. Were 
we to neglect mixing in other sectors, this expression would be valid by replacing 1/4 —>■ 
l/2C(rep). In the case of the wino-zino multiplet, which transforms as the adjoint, the 
appropriate replacement is 1/4 —> 1. But in addition one has to include the all important 
effect of mixing, and the corresponding formula is just as in eq. ( |6.11|) . This involves, in 
addition, 

/° = lim(QV^^)^4(mi,m2,g) 

q^Q 

— 4mim2/M^(m2 — m^)[m2 Inm^/ fJ.'^ — m^ Inm^///^ — {rn^ — nii)] (6.14) 

It should be kept in mind that the fermions in the neutralino sector are Majorana. 
The equations above apply to them too, provided one remembers to include factors of 1/2 
as appropriate. 

Explicit expressions for the integrals appearing in eqs. ( |6.1| ) and (6.8) are easily writ- 
ten. We refrain from doing so because the expressions are lengthy. A detail discussion of 
the behavior of these functions can be found in section 8 below. 



7. Top and Higgs 

The experimental extraction of gauge couplings in the standard model depends on 
assumptions on two other undetermined parameters, the top and Higgs masses. It is 
customary to make a global fit of observables including not just the gauge couplings but 
these masses — or at least the top mass — as well. How are we to deal with this in our 
translation between standard and extended models? 

Consider first the dependence on the Higgs mass. In supersymmetric models one 
must introduce two distinct Higgs doublets to ensure that all quarks and charged leptons 
are massive. Both fields get vacuum expectation values, and the physical spectrum is 
more complicated than that of the standard model. The details are spelled out in the 
next section. For now we would like to concentrate on the question of how to account 
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for the different Higgs sectors in the standard model and the supersymmetric one. The 
problem is that the physical Higgs of the standard model generally has no simple direct 
correspondence in the supersymmetric extension. The linear combination of fields that 
gets a vacuum expectation value, a clear candidate for the corresponding standard model 
Higgs, is not a mass eigenstate. 

Thus, in computing 5gk one must be careful to include these distinctions into AF^ in 
eq. (|3.4| ). In fact, AF^ should include all the supersymmetric Higgs sector contributions and 
subtract the standard model ones. However, the standard model calculation introduces an 
unknown parameter, the Higgs mass. In fact, some of the uncertainty in the determination 
of the standard model's gauge couplings is due to the lack of determination of the Higgs 
mass. In our approach, the standard Higgs contribution would be replaced by a sum of 
contributions from the two neutral scalars, plus, in addition, genuinely new contributions. 
When the two neutral scalars are degenerate and with mass equal to that of the standard 
model's Higgs, this sum is exactly equal to the standard model Higgs contribution that it 
replaces. Because the mass of the Higgs is unknown, and since the standard model gauge 
couplings are not very sensitive to it, it is a good approximation to retain the standard 
model Higgs contribution, and instead to neglect the corresponding terms in the model 
with two Higgs doublets. 

To summarize, given the uncertainties associated with the standard model determina- 
tion of gauge couplings, the calculation of Higgs sector threshold corrections needs include 
only contributions from physical scalars in the W and Z self-energy diagrams of fig. ^ 

It must be noted here that in the run-and-match analysis this is accomplished by 
neglecting mixing between the two Higgs doublets. Then a correspondence is made be- 
tween one of them and the standard model's, while the other gives a new light threshold 
correction. 

The contribution from the top quark is identical in the standard model and its exten- 
sion. It is included in the function Fi{gk) of eq. ( p.4|) . The problem is that the determina- 
tion of gauge couplings usually uses this as a fit variable. This could lead us into a com- 
plicated analysis of correlations. We are fortunate in that the models we are interested in 
— supergravity unification with radiative-induced electroweak symmetry breaking — yield 
top quark masses in the range obtained from these fits. Therefore, to good approximation, 
the effect of top need not be included in our calculations. 
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8. Review of the minimal SU{5) supergravity model. 

In the next section we wih ihustrate our treatment of the hght thresholds. In prepara- 
tion we review here the minimal SU{5) supergravity model. We also describe the procedure 
we have adopted for the calculation of the spectrum. The particle content of the model 
consist of three generation of quarks and leptons, the SU (3) x SU (2) xU{1)y gauge bosons, 
two Higgs doublets and the supersymmetric particles of all these particles. We neglect the 
Yukawa couplings of the two light generations. The superpotential of the effective low 
energy theory thus has the form 

W = \tU^Qh2 + XhD^Qhi + KE^Lhi + ^hih2 (8.1) 

where Q and L are quark and lepton doublet superfields, and U^, and are the corre- 
sponding SU{2) singlets; hi and /i2 are the two Higgs doublets. The soft supersymmetry 
breaking terms are 

^soft = iXtAtQU^h2 + XbAbQD'^hi + XrArLE^hi + B^ihih2 + h.c.)+ 

m^J/iil^ + m^J/i2|^ + m\\L\'^ + m\.\E''\'^ + mQ\Q\^ + mua\U''\^ + ml,.\D''\'^ . 

(8.2) 

The SU{5) model is specified by the boundary conditions at the unification scale. 
The low energy parameters can then be determined from knowledge of their scale depen- 
dence. This is dictated by the renormalization group. The renormalization group equations 
(RGEs) for the minimal supersymmetric standard model are well known and are given in 
appendix A for completeness. 

The boundary conditions at the unification scale that we adopt equate all of the 
trilinear coefficients Ax and equate the hard scalar masses in the soft supersymmetry 
breaking potential of eqn. ( |8.2| ). In addition one has conditions that follow from minimal- 
SU{5): equate the tau and beauty Yukawa couplings, Ab(MGUT) = At-(Mgut) and equate 
the gaugino masses Mi {i = 1, 2, 3, corresponding to the U (1) x SU (2) x SU (3) gauginos). 
The free parameters are then the common trilinear coefficients A, the soft Higgs mass 
parameter S, the top Yukawa At, the bilinear Higgs mixing coefficient |U, the common 
gaugino mass mi, the common scalar mass mo, the scale of unification Mqut and the 
unified gauge coupling acuT- The free parameter Ar is of course fixed by the tau mass. 

Given values of these free parameters one can compute the low energy spectrum of the 
theory, and then complete the calculation of the relation between gauge couplings at low 
energies as outlined in previous sections. We would like to explore the span of predictions 
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that correspond to a large region of this free parameter space. It is quite unnecessary, 
though, to vary freely over values of all of these parameters. In computing the spectrum, 
we take as initial values sin^ 6w = 0.233 and = 0.120, from which we determine Mqut 
and acuT- 

The Yukawa couplings of the heavy generation are given by 

1 nib T . ^ ''^t 

M,T = n — — M = - — n — (8-3) 

cos p V sm p V 

where v = \/v\-V v\ = 246/v^ GeV, tan/3 = v^jvx and v\^2 refer to the expectation 
values of the two Higgs doublets. We obtain the heavy generation Yukawa couplings by 
eq. (|8.3|) from the input parameters nit and tan/3, and the input values of m^- = 1.78 GeV 
and mf, = 4.0 GeV at the Z-scale. In our numerical analysis, we evolve the Yukawa 
couplings from the Z-scale to Mqut, and constrain the parameter space by the requirement 
Af,(MGUT) differs from Ar(MGUT) by no more than 5%. In so doing we neglect the effect 
of thresholds on the running of these Yukawa couplings. Thus for this procedure we need 
also a priori knowledge of the angle /3. 

In the next step in our numerical computation of the spectrum we evolve the Yukawa 
couplings, the soft SUSY breaking parameters and the gauge couplings with the boundary 
conditions specified at Mqut, to the electroweak scale, with the RGEs of the MSSM as 
given in appendix A. We then obtain the supersymmetric spectrum at the weak scale from 
the equations given below. 

In order to span over parameter space at the scale of grand-unification, and compute 
the low energy spectrum via the RGE's subject to the conditions that we obtain the correct 
masses of the tau-lepton and bottom-quark and that the correct electroweak breaking scale 
is the minimum of the tree level neutral Higgs potential, we trade the input parameters \\x\ 
and B for v and tan/3 (or, equivalently, v\ and v^)- Therefore, B and become computed 
parameters. This we can do because the running of the SUSY parameters, as given by the 
RGE's in appendix A, does not depend on the values of S or (except for \i itself). Thus, 
in our numerical analysis we span over the space generated by A, mt, mi , mo, tan/3 and 
the sign of /x, and compute the spectrum for each set of values. 

We turn now to a description of the spectrum itself. We lay down here the basis for 
our numerical treatment of the next section, making no attempt to describe the spectrum's 



salient qualitative features. A nice account of this can be found in ref. while extensive 



numerical analysis of the spectrum can be found in ref. fTj]. Our results for the spar- 
ticle spectrum agree qualitatively with the results of these references. For the two light 
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generation sparticles we neglect the Yukawa couplings in the RGE's. The light-generation 
sparticle masses may then be analytically calculated from the one-loop RGEs in terms of 
the three unknowns mi, mo, and cos 2/?: 

m| = mg + Cpm| + dp. (8.4) 

The coefficients cp for the different sparticles are 

Cp = csirrip) + C2(mp) -|- Y~Cy{mp) (8.5) 

with 

'=3(m,) = -?('l-(l + ^)"j (8.6a) 

where aouT is the coupling at the unification scale and t = log{mp/ Mqut). The RGE's 
are integrated to the physical sparticle mass to obtain the c's. The d's arise from D-terms 
in the potential and are given by 



dp = 2{ - tan^ 9w ] cos 2/3 . (8.7) 



In case of small tan/3 the bottom quark and tau lepton Yukawa couplings can be 
neglected as well. In the case of large tan /3 the heavy generation sparticle mass matrices 

arc given by 



^ I ml + m'^ + mt{At + cot /?) \ 

^ \mt{At + ^ coi fi) mfj, + + d^ 

M?=( ^Q + '^b+dd, mb{At + n tan /3)\ 

^ \mb{Ab + n tan (3) fnjyc + ml + d^^ J 

M? = ( ^l+^T + dei mr{Ar + ^itanP)\ g , 

\mr{Ar + IJ, tan P) m\c + ml + d^^ J \ ■ ) 

The mass eigenstates are obtained by diagonalizing the mass matrices above by a 2 x 2 
unitary transformation. 
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The chargino and neutralino mass matrices are given by 

M2 MwV^ sin P 



and 



/Ml ~Mzsm9wcosP Mzsm9wsmP \ 

M2 -^z cos 9w sin — cos 9w cos /3 

— sin 9w cos /3 cos 9w sin /3 ' 

V sin 6*^1/ sin /3 —M^ cos 6*^1/ cos /3 / 

(8.10) 

respectively, where Mi = {ai/ajj)mi (i = 1,2,3) are the gaugino masses at their mass 
scale. The mass eigenstates and mixing can be found in ref. [^. 

There are two complex Higgs doublets in the MSSM, (/ij, /i2), whose neutral vacuum 
expectation value (VEV) V2 gives mass to the | charged quarks, and (/i^, h^) whose neutral 
VEV gives mass to the charged quarks. Of the eight real degrees of freedom, three are 
the Goldstone bosons that are absorbed by the standard model gauge bosons. The CP- 
odd neutral Higgs boson and the orthogonal combinations h'^ = sin (3(h^)* + cos/3 /ij, 
h~ = (h^)* correspond to physical particles with masses 

= -m3(tan/3 + cot/3) , (8.11) 

and 

ml± = ni^ + m\ , (8.12) 

where ma = Bfi. The two CP-even neutral Higgs bosons mass eigenstates are obtained 
by diagonalizing 

2 _ sin2^ /cot /3M| + tan /3m^ -(M| + m^) \ , . 

2 V -(M^z+^a) tan/^Ml + cot/^m^+Ay) ' ^^-^^^ 

where A is the one loop radiative correction to the neutral Higgs boson masses 

' log + (8.14) 



a 9 2 o I 2 

87r^ \ mf 

The linear combination of Higgs fields that gets a VEV, hsrm and the one that does not 
get a VEV, hgusy, are obtained from the mass eigenstates h and if by a rotation, 

hsm \ ^ [ cos(a - (3) sin(a - P) \ f h\ , . 

hsusyj \-sin{a-P) cos{a~P) J \H J ' 
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where a is the angle that rotates from the gauge to the mass eigenstates. Requiring a 
negative eigenvalue for the neutral Higgs mass-squared matrix and that the Higgs potential 
is bounded from below imposes two conditions on the running of the mass parameters: 

1. m\m^ — mg < 

2. + 7712 - 2|m3| > 

where 2 = 2'^ ^ "^3 ~ 

9. Numerical Results and Discussion 

We first consider the functions $x of section 5 as compared with the naive (leading- 
log, or "run-and-match") result. The latter can be obtained by the standard run-and- 
match computation, or simply by taking the large m/Mz limit of our expressions. Both 
procedures must agree when in our expressions SU{2) x U{1) breaking is neglected. 

Consider the contributions of scalars to ^ught of eqs. ( |5.11| ) and (|5.13| ), both in our 
procedure and using the run-and-match approximation. When SU{2) x U{1) breaking 
is neglected, so that all members of a weak-isospin multiplet are degenerate, all three 
functions $x are simply proportional to the fimction Ii . Thus, for the contribution of any 
one scalar, the ratio of 5iight computed in the run-and-match approximation to our result 
is simply the ratio of the leading- log term in Ii to the whole expression for Ii. This ratio is 
plotted in fig. ^ as a function of m/Mz- It is apparent that the error in the approximation 
is small except for m/Mz ~ 1- 

Next, consider the contributions of fermions to ^ught- The $g piece is easiest to 
analyze, since it is simply proportional to the function I2; see eqs. ( |6.8q| ) and ( |6.9|) . In 
fig. P we plot the ratio of the naive run-and-match approximation to the complete one 
loop threshold correction function $3, due to virtual fermions. In this case, again, the 
convergence to the naive run-and-match limit is extremely fast. 

The most interesting case is that of the fermion contributions to $z and ^w- Recall, 
from section 6 that $z and $VK are linear combinations of the integrals I2, I3 and 14. The 
behavior of I2 was analyzed above, and was found to be quite uninteresting. Turning our 
attention to /s, we see by inspection of eq. ( |6.8£| ) that its asymptotic behavior, for large m, 
is /3(m, m, q) ~ —Am?/q^ ln(m^/ /U^) + 2/3 + /Ibm? -|- O^q'^ /m^). This is an example of 
non- decoupling, in which the leading term grows as a power (rather than a log) of the mass. 
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In fact, the combination — $VK is precisely what enters the computation of Veltman's p 
parameter, which is weU known to grow as ~ ml in the standard modelll^. Notice though 
that for particles in vectorial representations the prefactor (r£ — r^)^ vanishes. So a large 
effect will depend on the breaking of SU{2) x U{1). 

It is interesting to note that the combination /s + /4 which often enters into $2, 
vanishes in the limit of degenerate multiplets. This can be seen by writing it as 

13(1711, 1712, q) + 14(^1,1712, q) = — 0(1712 -mi) / dx {xm2 - (l - x)mi) log M'^ (x) / . 

r Jo 

(9.1) 

These type of leading corrections to 5iight are missed altogether by the naive run- 
and-match calculation. In fact, there is no pure log term in I3, an indication that it is 
invisible to the naive approach. And, clearly, the combination I3 + is missed altogether 
too. Therefore we cannot quantify the effects of these terms, relative the run-and-match 
approach, in a manner similar to the above examples. We postpone further consideration 
until we come to computations in the full SU{5) supergravity model, below. 

To get some idea of how important the SU{2) breaking effects can be, we can look at 
the contribution to ^ught from a doublet of scalars, with masses mi and m2 respectively, 
and compare to the would be contribution from a degenerate doublet with some average 
mass, mavg- For definiteness we shall take m^vg = {rni + m2)/2, although one could just 
as well use, say, m^vg = ^/rr^m2. In fig. ^we plot the ratio /i(mi, m2, q) / Ii{rnavg, rriavg: q) 
for m2/q = 1.5, 3.5, 5.5 and 7.5 (corresponding to the ever increasing location of the peaks 
in the graph). It is clear that SU{2) breaking effect can be rather considerable. It should 
be pointed out that large SU (2) breaking effects would show up as significant deviations 
from unity of Veltman's p-parameter. For consistency with the observed value of p, the 
parameters of the SUSY-GUT are such that the SU{2) breaking effects are small. 

These comparisons, although instructive, are incomplete. They miss the possibility of 
cumulative effects of the many particles in the rich supersymmetric spectrum. We therefore 
use our equations in the context of the model described in the previous section. This will 
allow consideration of the effects of mixing. These are obviously important and altogether 
missed by the run-and-match approach, but serious consideration of them requires a choice 
of reasonable mixing. In the minimal supergravity model based on SU{5) the mixing is 
completely fixed in terms of other more fundamental parameters. 

We calculate the spectrum in the minimal supersymmetric SU{5) model as described 
in section 7. There are five input parameters, namely A, mo, mi, tan/3 and the top 
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quark mass rrit. In addition one must specify the sign of n. They are constrained by the 
requirement of radiatively induced symmetry breaking of the electroweak group, and lower 
bounds on supersymmetric particle masses. We perform a systematic, although rather 
coarse-grained, search through parameter space. Each of the five parameters is allowed to 
vary over a reasonable range. The range is sampled at some fixed predetermined interval 
size. Table 1 summarizes the range and size of the sampling interval for each of the 
five parameters. This represents a sampling of more than 4 x 10^ points. Of the resulting 
spectrum we require that it does not violate experimental bounds: we require that chargino, 
squark, gluino, slepton and neutralino masses to be in excess of 45 GeV, 100 GeV, 150 
GeV, 43 GeV and 20 GeV, respectively. We also require that the lightest neutralino is 
lighter than the chargino. Only some 3 x 10* sampling points in parameter space survive 
this requirements. To study the dependence on each of the input parameters we use a 
finer grid for one of the inputs and take two or three extreme points for the other input 
parameters. We find that rrit < 160 GeV and tan/3 < 30. This is the case in part because 
we have assumed that the tau and bottom Yukawa couplings are equal at the unification 
scale. The gluino mass limit, Mg > 150 GeV, constrains mi to be above 55 GeV. Upper 
limits on \Ao\, uiq and mi can be set by resorting to naturalness arguments 0. However, 
as our interest is in the SUSY spectrum close to the electroweak scale, we impose upper 
bounds on these parameters not far above the electroweak scale. 







Xf 


AX 


Ao 


-200 


200 


50 


mo 





400 


50 


mi 

2 


100 


300 


25 


tan/3 


2 


30 


4 


mt 


110 


170 


10 



Table 1. The range and sampling size of the parameter space. Each free parameter X is 
sampled in the interval {Xi.Xf) with spacing /S.X between consecutive points. 
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How good is the run-and-match analysis in this case? In fig. |^ we compare our 
calculations, in (a) method-I, and (b) method-II, with the result of ref. 0. We plot sin^ 6*, 
calculated in (|2.5| ) and (|5.11| ) or (|5.13| ) using as = 0.120 and Oe.m. = 1/127.9, against sin^ 9 
in the run-and-match approximation. Each point on this plot represents a specific choice 
of the free parameters. It is clear that for many points the run-and-match approximation 
is far too crude to be useful. The effect of two-loop running can be incorporated in these 
graphs (and the following ones) by a common shift of all points. The two band structure 
that can be seen in fig. || come from the two signs of fj,, the lower higher values corresponding 
to positive fi. One should not worry much about the discrepancy between the two graphs 
in fig. § since they correspond to different interpretations of the same "experiment" . 

In figs. p|-p!2| we plot the dependence of sin^ 6 in our calculation on the mass of a few 
of the superparticles, in (a) method-I, and (b) method-II. The dependence on the lightest 
chargino mass, min, can be seen in fig. ^. As expected, the points tend to concentrate 
around the sin^ 9 value of the run-and-match approximation as Mcmin increases. The mass 
of the lightest neutral Higgs increases rapidly as a function of tan/3, reaching a limiting 
value of about 100 GeV at tan/3 ~ 6. This, rather than the two signs of /i, is the reason 
for the two band structure of fig. ^ a scatter plot of sin^ 6 against Mh, min- The band 
at about 60 GeV is made up mostly of points at tan/3 = 2, while the band at 100 GeV 
consists mostly of points with tan/3 > 6. Fig. 11 shows sin^ ^ against the mass of the 
lighter of the two top-squarks. It is similar to fig. ^ in that points tend to accumulate at 
the run-and-match value as MT,min increases, although the spread seems larger. Fig. 12 
again demonstrates these effects, now as a function of the lightest neutralino. Finally, 
in fig. |13| we verify that 5p remains small in method-I. Note that the larger values of 5p 
correspond to a narrow band of sin^ 9. 

We would like to emphasize that our numerical analysis is not intended as a complete 
analysis of the parameter space and of the SUSY spectrum. Nor do we believe, in view of 
the huge parameter space, that such analysis carry substantial meaning. The parameter 
space can be constrained by making specific assumptions, motivated to some extent by 
theoretical considerations, on the form of the Kahler potential in the supergravity model, 
and consequently on the initial boundary conditions at the unification scale. Partial anal- 
ysis of the SUSY spectrum, based on such specific assumptions, have been performed by 



several groups[|IJ]. In this paper, our purpose is to address the question of how to appropri- 
ately account for the effects of light thresholds, in a way that can consistently incorporate 
two-loop running, and to illustrate the need for a proper method to do so. Our numerical 
analysis should be regarded merely as illustrative. 
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10. Conclusions 



The standard run-and-match approach to the calculation of light threshold effects 
is seen to be inaccurate when the supersymmetric spectrum is not too heavy and there 
is substantial SU{2) breaking and mixing. We have described a better approximation 
which incorporates these effects. Moreover, our method incorporates readily the effects of 
two-loop running. Although we have concentrated on the predicted value of sin^ 9 in the 
numerical analysis in this work, clearly similar analysis apply to the calculation of other 
quantities of interest. For example, one may include threshold effects in the calculation of 
the spectrum itself [jl^ . 

The calculation that we have presented is a sensible approach, as long as supersym- 
metric particles are not found and provided there is no evidence of experimental deviations 
from standard model expectations. When either of these assumptions is violated the cal- 
culation should be replaced by a full calculation of the novel process, including possible 
non-universal effects, that is, non-oblique corrections. 

In the mean time one may use this type of calculation to restrain more severely the 
allowed parameter space of supersymmetric grand-unified theories. It would be interesting 
to produce an analysis including simultaneously the predictions of sin^ 9 and other re- 
stricting quantities, to better limit the boundaries of parameters space. For example, the 
inclusive rate for radiative B decays may be fairly restrictive of extreme values of tan/3. 
Unfortunately, to the extent that the standard model is successful, one can always find 
a region of the SUSY-GUT that satisfies all experimental constraints, for the sparticles 
decouple as they become much heavier than the massive vector bosons, and in that limit 
the predicted value of sin^ 9 is consistent with experiment. 
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Appendix A. RGE's for the MSSM 

The renormalization group equation for the minimal supersymmetric standard model 
are well known and are given here for completeness. 



(3A? + \\l - - 
(3A^ + ^A? + - 

(2A^ + Ixl - \gl - 

{2,\lAb + \\^,A, + \\IA, + ^-glM^ + ^glM^ + ^^/i^i), (A-l/) 

{2\lAr + \\lA^ + \glM2 + ^^?Mi), (A.l^) 

(^A?A, + ^A^A + \xlAr + ^^1^2 + ^glMi). (A.l/^) 

(3n + i^r - 3^7iM| - \glMl), (A.li) 
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16 ~ _ 1 _ 
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(2F. - fgfM!), (A.lo) 
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where, 

Ft = Xl {ml + m^c + ml^ + A^) (A.2a) 

Ffe = Xliml + ml. + m^^ + Al) {A.2b) 

Fr = Xliml + ml. + ml^ + Al) (A.2c) 
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Figure Captions 

Fig. 1. Tree level Feynman diagram for the scattering amplitude for |Uz/^ eue scatter- 
ing. 

Fig. 2. Feynman diagrams for one loop virtual effects of new heavy particles. 

Fig. 3. One-loop Feynman diagrams in the calculation of threshold corrections: scalar 
particle contributions. 

Fig. 4. One-loop Feynman diagram in the calculation of threshold corrections: fermion 
contributions. 

Fig. 5. Ratio of the naive run-and- match approximation to the complete one loop thresh- 
old correction functions due to a virtual scalar. 

Fig. 6. Ratio of the naive run-and-match approximation to the complete one loop thresh- 
old correction function $g, due to virtual fermions. 

Fig. 7. The effect of SU{2) breaking, as included in our analysis, in the scalar contri- 
bution to ^w- This is a plot of the ratio /i(mi, m2, q) / Iiijn.^^^, ^i^avg, q), where 
"^avg = ("^1 + "^2)/2, for m2/q = 1.5, 3.5, 5.5 and 7.5 (corresponding to the ever 
increasing location of the peaks in the graph) . 

Fig. 8. Plot of the calculated sin^ 9, this work, in (a) method-I, and (b) method-II, ver- 
sus the calculated sin^ 6* in the standard 'run-and-match' approximation 0. See 
eq. ( p.5|) . Each point on this plot represents a specific choice of free parameters 
of the minimal super symmetric SU{5) grand-unified theory that survives exper- 
imental bounds on the spectrum. The parameter space spanned here is detailed 
in Table 1. 

Fig. 9. Scatter plot of sin^ 6 against the mass of the lightest chargino, Mcmim for the 
minimal supersymmetric SU{5) grand-unified theory, in (a) method-I, and (b) 
method-II. Each point represents a choice of parameters; see Table 1. 

Fig. 10. Scatter plot of sin^ 6 against the mass of the lightest neutral Higgs, MH,min, for 
the minimal supersymmetric SU{5) grand-unified theory, in (a) method-I, and 
(b) method-II. Each point represents a choice of parameters; see Table 1. 

Fig. 11. Scatter plot of sin^ 9 against the mass of the lightest top-squark, MT,min, for the 
minimal supersymmetric SU{5) grand-unified theory, in (a) method-I, and (b) 
method-II. Each point represents a choice of parameters; see Table 1. 

Fig. 12. Scatter plot of sin^ 9 against the mass of the lightest neutralino, M^r min, for the 
minimal supersymmetric SU{5) grand-unified theory, in (a) method-I, and (b) 
method-II. Each point represents a choice of parameters; see Table 1. 
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Fig. 13. Scatter plot of sin^ 9 against the deviation of Veltman's p parameter from unity, 
for the minimal supersymmetric SU{5) grand-unified theory in method-I. Each 
point represents a choice of parameters; see Table 1. 
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